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Abstract. We have studied isospin breaking effects in the X(3872) resonance and found a natural explanation for the 
branching fraction of the X decaying to // iff with two and three pions being close to unit. Within our framework the X(3872) is 
a dynamically generated resonance in coupled channels. We also study the relationship between the couplings of the resonance 
to the coupled channels with its wave function, which further helps us to understand the isospin structure of the resonance. 
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& 1. INTRODUCTION 



The X(3872) resonance has been first observed by the Belle collaboration [1] and later confirmed by other four 
experiments the CDFII, DO and BaBar collaborations [2, 3, 4]. It has first been seen in the decay channel // y/n + 7t~ and 
the analysis of the nn spectrum indicates that this meson pair comes from a p meson. There has been no observation 
of a charge partner for the X(3872) state, which indicates an isoscalar nature for it. Later on the X(3872) has been 
observed in another two decay channels [5]: The channel J/y/y tnat fi xes me C-parity of the state as positive and 
the channel J/ \]/n + n~n° where the three pion spectrum shows a peak for an CO resonance. The measured branching 
fraction 



Si ^Ir J 'T +i r^ = 1.0±0.4 ± 0.3, (1) 

raises a problem, since it seems to indicate a huge isospin violation in the decays of the X(3872), since the case of p 
production (two pions) has isospin 1 while the CO has isospin 0. A careful analysis of the data on this state has been 
5— 1 ' done in [6], concluding that its quantum numbers must correspond to J p — 1 ++ or J p = 2 h . 

From the theoretical point of view the X(3872) mass is very close to the D°D*° threshold and it has been suggested 
[7, 8, 9, 10] that this state is a DD* molecule. There has been though some discussion on whether the D + D* 
component is important in the structure of the state, since this component is bound by about 8 MeV, while the neutral 
component is bound by less than 1 MeV [11, 12]. 

In this work we explain how the branching fraction of Eq. (1) is naturally explained for a dominant isospin nature 
of the X(3872). We also stress the importance of taking into account the charged DD* components in the structure of 
the X even though the neutral wave function extends much further in space than the charged one. This is a consequence 
of the fact that only the values of the wave function around the origin are important in the description of short range 
processes as the decays to J / y with two or three pions. 



2. THE DECAYS OF THE X(3872) INTO TWO AND THREE PIONS 

We follow here the approach of [12] where the X(3872) is generated dynamically within coupled channels. As noted 
in [10, 12] only the D + D* and D°D*° components of the state provide sizeable couplings and therefore we consider 
only these two channels in the approach. 



The model is based on the calculation of the unitarized transition matrix (T-matrix). For the tree level interaction 
of the mesons, Lagrangians based on the SI/ (4) extensions of the hidden gauge Lagrangians are used. The SI/ (4) 
symmetry is broken in order to take into account the bigger mass of the charmed and hidden charmed mesons. The 
tree level amplitudes are used as kernel to solve a scattering equation that in an on-shell approach takes the form 



T = {l-VG)~ l V (2) 

where T is the T-matrix, V is a matrix containing the tree level amplitudes projected in s-wave and G is a diagonal 
matrix with each diagonal element containing the two particle propagators for each channel. 

The loop functions in G are divergent and can be regularized using dimensional regularization, which gives a 
subtraction constant a as a free parameter that can be used to adjust the pole generated in the T-matrix of Eq. (2) 
to its right position, fitting in this way the X(3872). 

Close to the pole the T-matrix can be written in the separable form 



Tij = 8 ^ 8i (3) 

s Spole 

where s is the invariant mass squared, s po \ e is the pole position and are the couplings to each channel. Using Eq. (3) 
one can calculate the couplings of the resonance to each channel by calculating the residues of the T-matrix at the 
pole position. By using physical masses for the charged and neutral components of the D-meson doublet members, 
one finds a difference of at most 1.5% in the couplings of the X to the charged and neutral channels. This is a tiny 
difference which we will neglect in what follows. 

In order to calculate the ratio of the decays into J / y/p and J/y/co one has to take into account the coupling first of 
the X to a D°D*° or D + D* and then the coupling of the DD* components to J/ij/p or J /yco. In the hidden gauge 
formalism this former coupling is equal for both cases. What differentiates the situation with the p meson from the 
situation with the CO mesons is that in the case of the p which has isospin 1 the loops with neutral and charged DD* 
mesons interfere destructively while in the case of the CO with isospin the interference is constructive. As a result the 
ratio of the coupling squared of the X to p and CO is given by 



Rp/co = ^-^=0.032 (4) 

Gll+Cr22 

The value 0.032 was obtained in the limit of binding energy with the loops calculated with dimensional regular- 
ization. 

In order to calculate the full ratio of the branching fractions one still has to take into account the phase space 
available for the p to decay into two pions and the co into three pions. For that we integrate the spectral function of 
each meson taking into account the much bigger width of the p and the kinematical constrains. The result obtained is 



@(X^J/n,K+K-K°) _ 

@{X -^J/yn+n-) 

which is compatible with the value 1 .0 ± 0.4 from experiment [5] . One sees that although the couplings of the X to the 
isospin 1 state J /y/p is suppressed in relation to the isospin state J/xj/co by a factor 30, the much bigger phase space 
for the p to decay than for the co, compensates this suppression and brings the branching fraction to a value close to 
the observed experimentally. 



3. THE WAVE FUNCTIONS OF THE X(3872) 

To calculate the wave functions of the X resonance we start from a phenomenological potential that in momentum 
space reads: 



(p'\V\p)= V(p',p)= v&(A-p)&(A-p') 



(6) 



where the cut off A is fixed in order to find the resonance with the appropriate mass. To do that we write the T-matrix 
and the Lippmann-Schwinger equation: 



(p\T\p') = T(p',p)=t&(A-p)&(A-p>) 
t = (l-vG) _1 v 



the coupled channels loop function G here is a diagonal matrix given by: 



G = 



Gn 

G 22 



, G tt = [ 



d 3 p 



P<^E-M: 



(7) 
(8) 



(9) 



The poles of the T-matrix are given by the equation det(l — vG) = 0. The value of Lambda is chosen so that this 
equation is satisfied for the loop calculated at the energy where we want to have the X state bound. 
We write then the Schorodinger equation for the potential: 



(H + V)\w) = E a W) 



(10) 



where Hq is the free Hamiltonian and E a is the energy of the resonance. Projecting this equation in momentum space 
one obtains a set of two coupled equation for the wave function in each channel: 



&(A-p) 



(<%!> + <*|lfc>) 



(p\\l/2) = v- 



®(A-p) 



w f d\({k\ wl )+{k\ w2 )) 

E— Jk<A V / 



(11) 
(12) 



The two integrals in the right hand side of the equations are constants that can be calculated by normalizing the wave 
function. 

The couplings of the state to the channels can again be calculated from the residues of the T-matrix: 



2 2 _ 2 
81=82^8 



Jim (E-E a )tij = - 



dGn dG 2 2 \ 1 
dE dE J 



(13) 



E=E a 



In the limit when E^ ->■ 0, we have ^ 



E=E a 



2 2 
8l=82 



00 and we find 



7i 



-Bl 



where 



L Bi 



Mi-E a . 



(14) 



(15) 
(16) 



and the M, is the value of the threshold for channel i. 

As shown in [13] working out Eq. (11), Eq. (12) and Eq. (13) one obtains 



*<??!=/ d 3 p(p\ Vl ) 

J p<A 

gG% 2 = [ d 3 P (p\y 2 ) 



(17) 
(18) 



TABLE 1. Comparative results for different potentials for a D°D*° binding energy 
of 0.1 Me V. 



Form 


A 


FT 

g 


Vl(0)/V2(0) 


W1/Y2 




V>2 


R p/a> 


Factor 


[MeV] 


[MeV] 













Sharp 


653 


| 3202 [ 


1.31 


I 1.31 


| 3.29 


2.50 


0.018 


Gauss 


731 


j 3238 ( 


1.20 


I 1.29 


| 3.30 


2.56 


0.016 



Lorentz I 834 I 3254 I 1.17 I 1.28 I 3.31 I 2.58 I 0.015 



and we can recognize that the integrals in these equations are the Fourier transform of the momentum wave function 
for r = 0. So we can rewrite Eq. (17) and Eq. (18) as 

gG? 1= (2703/2^(0)= ^ (19) 

gG% 2 = (27T) 3 /V 2 (6) = v> 2 (20) 

These equations show that the couplings of the the resonance to its channels are proportional to the value of the 
wave function at the origin (denoted as i/a). 

Had we chosen another kind of form factor for regularize the potential in Eq. (6) we would arrive at the same 
equations, but with a slightly different definition of \jf: 

Wi = J d 3 kf(k)(k\ Wi ) (21) 

where f(k) is a form factor substituting the ©(A — k) in Eq. (6). This new form factor has also a scale A that in 
coordinate space represents a range about 1/A. Since a typical cut off of 700 MeV represents in coordinate space a 
range of about 0.3 fm, y represents here a smeared value of the wave function around the origin. 

In Table 1 we show results for potentials with three different form factors. The wave functions for each kind of 
potential are slightly different but the values of i//, and the ratio of these values do not change much, since these are 
the quantities that carry the information about the short range physics. 

In Figure 1 we show plots of the wave functions and probability densities for the sharp cut off potential and plots 
comparing the wave functions for the three different kinds of form factors. 

With regard to the isospin structure of the X (3872) we can first write its wave function in terms of the charge basis 
and then transform it to isospin basis: 



|X) = »,,(?)|D D*°) + «, 2 (?)|D+D-) 

We see from Eq. (23) that the mixing between the two possible isospin states depends on the relative distance of the 
two D mesons, r. But in physical processes this mixing will be given by the values of the wave function around the 
origin. Using the values given in Table 1 one sees that for short range processes the contribution of the isospin 1 
state is about 2% and the X can be regarded as an almost pure isospin state. 



(22) 
(23) 



4. CONCLUSIONS 

We have studied the isospin structure of the X(3872). In our model this state is generated dynamically in coupled 
channels and is interpreted as an s-wave bound state of DD* . The apparent isospin violation in the branching fraction 
ratio of X — > J/ and X — > J/\i/n + 7t~7t° is naturally explained due to the much larger phase space for the p to 
decay into n + n~ than the a> to decay into n + n^n . 

The mass of the X(3872) is much closer to the D°D*° threshold than to the D + D* which is bound by around 
8 MeV. As a consequence the wave function of the X in the D°D*° channel extends much further on space than in 




FIGURE 1. Neutral and charged Wave function components for a D°D*° binding energy of 0. 1 MeV. In the upper left panel 
one can see the value of the wave functions at the origin for both channels, in the lower left panel we plot the probability density 
for each channel. In the right panel we show plots of the wave functions for each channel calculated with the three different form 
factors in the potential. 

the D + D* channel. Nevertheless physical short range processes are dominated by the couplings of the state to each 
channel and these couplings have a negligible isospin violation. In terms of wave functions this can be understood 
since the couplings reflect an averaged value of the wave function around the origin and not in the whole space. 
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